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In this paper, we review our novel information geometrodynamical approach to chaos (IGAC) on 
curved statistical manifolds and we emphasize the usefulness of our information-geometrodynamical 
entropy (IGE) as an indicator of chaoticity in a simple application. Furthermore, knowing that 
integrable and chaotic quantum antiferromagnetic Ising chains are characterized by asymptotic 
logarithmic and linear growths of their operator space entanglement entropies, respectively, we 
apply our IGAC to present an alternative characterization of such systems. Remarkably, we show 
that in the former case the IGE exhibits asymptotic logarithmic growth while in the latter case the 
' IGE exhibits asymptotic linear growth. 

\ At this stage of its development, IGAC remains an ambitious unifying information-geometric 

theoretical construct for the study of chaotic dynamics with several unsolved problems. However, 
^ ■ based on our recent findings, we believe it could provide an interesting, innovative and potentially 

powerful way to study and understand the very important and challenging problems of classical and 
quantum chaos. 
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INTRODUCTION 



In classical and quantum dynamics there is no unified characterization of chaos. In the Ricmannian [1] and Finslerian 
(a Finsler metric is obtained from a Riemannian metric by relaxing the requirement that the metric be quadratic on 
each tangent space) geometrodynamical approach to chaos in classical Hamiltonian systems, an active field of research 
concerns the possibility of finding a rigorous relation among the sectional curvature, the Lyapunov exponents, and 
the Kolmogorov-Sinai dynamical entropy (i.e. the sum of positive Lyapunov exponents) ||. The largest Lyapunov 
I exponent characterizes the degree of chaoticity of a dynamical system and, if positive, it measures the mean instability 
t^-j- . rate of nearby trajectories averaged along a sufficiently long reference trajectory. Moreover, it is known that classical 
chaotic systems are distinguished by their exponential sensitivity to initial conditions and that the absence of this 
property in quantum systems has lead to a number of different criteria being proposed for quantum chaos. Exponential 
decay of fidelity, hypersensitivity to perturbation and the Zurek-Paz quantum chaos criterion of linear von Neumann's 
entropy growth Q are some examples [f|. These criteria accurately predict chaos in the classical limit, but it is not 
r* ■ clear that they behave the same far from the classical realm. 
. £h ! The present work makes use of the so-called Entropic Dynamics (ED) |6| . ED is a theoretical framework that arises 
from the combination of inductive inference (Maximum relative Entropy Methods, 7]) and Information Geometry 
(Riemannian geometry applied to probability theory) (IG) Q. As such, ED is constructed on statistical manifolds. 
It is developed to investigate the possibility that laws of physics - either classical or quantum - might reflect laws of 
inference rather than laws of nature. 

This article is a follow up of a series of the authors works @, El El El El El EE Eft Especially the work 
presented in [TEj will be discussed in more detail. The ED theoretical framework is used to explore the possibility 
of constructing a unified characterization of classical and quantum chaos. The general formalism of the IGAC is 
presented by investigating a system with 31 degrees of freedom (microstates), each one described by two pieces of 
relevant information, its mean expected value and its variance (Gaussian statistical macrostates) . This leads to 
consider an ED model on a non-maximally symmetric 6Z-dimcnsional statistical manifold A4 S . It is shown that M s 
possesses a constant negative Ricci curvature that is proportional to the number of degrees of freedom of the system, 
TZm s = SI. It is shown that the system explores statistical volume elements on M s at an exponential rate. We 
define an information geometrodynamical entropy (IGE) Sm s of the system and we show it increases linearly in time 
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(statistical evolution parameter) and is moreover proportional to the number of degrees of freedom of the system. 
The geodesies on M s are hyperbolic trajectories. Using the Jacobi-Levi-Civita (JLC) equation for geodesic spread, 
it is shown that the Jacobi vector field intensity Jm s diverges exponentially and is proportional to the number of 
degrees of freedom of the system. Thus, TZm 3 > $M 3 and Jm, are proportional to the number of Gaussian-distributed 
microstates of the system. This proportionality leads to conclude there is a substantial link among these information- 
geometric indicators of chaoticity. We emphasize that our IGE provides an information-geometric analog of the Zurek- 
Paz quantum chaos criterion [10| . As a physical application of our general theoretical scheme that we have called 
the Information Geometrodynamical Approach to Chaos (IGAC), we provide an information-geometric analogue of 
quantum energy level statistics for integrable and chaotic quantum spin chains. It is known [ItJ that in the integrable 
case, the antiferromagnetic Ising chain is immersed in a transverse homogeneous magnetic field and the level spacing 
distribution of its spectrum is the Poisson distribution. Instead, in the chaotic case, the antiferromagnetic Ising chain 
is immersed in a tilted homogeneous magnetic field and the level spacing distribution of its Hamiltonian spectrum 
is the Wigner-Dyson distribution. The antiferromagnetic Ising spin chain in external magnetic field is one example 
of order-to-chaos transition in quantum many-body context and it is used here as a demonstrating example of the 
conjectured connection between the Wigner-Dyson (Poisson) statistics and nonitegrability (integrability) in quantum 
mechanics. Moreover, it is known that integrable and chaotic quantum antiferromagnetic Ising chains are characterized 
by asymptotic logarithmic and linear growths of their operator space entanglement entropies [l7| , respectively. 

Following the results provided by Prosen, we study the information-geometrodynamics of a Poisson distribution 
coupled to an Exponential bath (regular case) and that of a Wigner-Dyson distribution coupled to a Gaussian bath 
(chaotic case). Remarkably, we show that in the former case the IGE exhibits asymptotic logarithmic growth while 
in the latter case it exhibits asymptotic linear growth. 

The layout of this paper is as follows. In Section II, the general formalism of the IGAC is applied to study a 
simple example, an ED Gaussian statistical model. In Section III, the main indicators of chaoticity within our novel 
theoretical construct are introduced by studying the ED Gaussian model. In Section IV, special focus is devoted to the 
role of the IGE as an indicator of temporal complexity on curved statistical manifolds. In Section V, after presenting 
the basics of the IG of Poisson and Wigner-Dyson distributions, we briefly review the conventional approach suitable to 
study the energy level statistics of integrable and chaotic antiferromagnetic Ising chains immersed in external magnetic 
fields. In Section VI, we present our IGAC-based novel characterization of the quantum energy level statistics of such 
Ising chains. Finally, in Section VII, we present our final remarks. 

II. THEORETICAL STRUCTURE OF THE IGAC: A SIMPLE EXAMPLE 

The IGAC arises as a theoretical framework to study chaos in informational geodesic flows describing physical, 
biological or chemical systems. A geodesic on a curved statistical manifold represents the maximum probability 
path a complex dynamical system explores in its evolution between the initial and the final macrostates. Each 
point of the geodesic is parametrized by the macroscopic dynamical variables defining the macrostate of the system. 
Furthermore, each macrostate is in a one-to-one relation with the probability distribution representing the maximally 
probable description of the system being considered. The set of macrostates forms the parameter space while the set 
of probability distributions forms the statistical manifold. The parameter space is homeomorphic to the statistical 
manifold. IGAC is the information-geometric analogue of conventional geometrodynamical approaches [3, 0] where 
the classical configuration space Te is being replaced by a statistical manifold A4$ with the additional possibility 
of considering chaotic dynamics arising from non conformally flat metrics (the Jacobi metric is always conformally 
flat, instead). It is an information-geometric extension of the Jacobi geometrodynamics (the geometrization of a 
Hamiltonian system by transforming it to a geodesic flow [HI]). The reformulation of dynamics in terms of a geodesic 
problem allows the application of a wide range of well-known geometrical techniques in the investigation of the solution 
space and properties of the equation of motion. The power of the Jacobi reformulation is that all of the dynamical 
information is collected into a single geometric object in which all the available manifest symmetries are retained- the 
manifold on which geodesic flow is induced. 

Using information-geometric methods, we have investigated in some detail the still open problem of finding a 
unifying description of classical and quantum chaos [ic| |. One of our goals in this paper is that of representing an 
additional step forward in that research direction. 

A. The ED Gaussian Model 

Maximum relative Entropy (ME) methods are used to construct an ED model that follows from an assumption about 
what information is relevant to predict the evolution of the system. Given a known initial macrostate (probability 
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distribution) and that the system evolves to a final known macrostate, the possible trajectories of the system are 
examined. A notion of distance between two probability distributions is provided by IG. As shown in [Tg . [20| this 
distance is quantified by the Fisher-Rao information metric tensor. 

In the following example, we consider an ED model whose microstates span a 3Z-dimensional space labelled by the 



) 2 > x 3 



a = 1, 



Z and 



» 



with 



1, 2, 3. We 



variables |x| = f( 2 \...., z«} with f( Q ) = (x[ a) 

assume the only testable information pertaining to the 31 degrees of freedom j consists of the expectation values 
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and variances A,<«> s J((4°> - (4"'» 2 ). The set „, these expectation valnes define the 6,-d«io„al 

space of macrostates of the system. A measure of distinguishability among the macrostates of the ED model is 
obtained by assigning a probability distribution P (^X\Q^j to each macrostate where {©} = {^'^i"': ( ' 2 '^i Q ' ) | 
with a — 1, 2,...., I and a = 1, 2, 3. The process of assigning a probability distribution to each state endows A^s 
with a metric structure. Specifically, the Fisher-Rao information metric is a measure of distinguishability among 
macrostates. It assigns an IG to the space of states. Each macrostate may be viewed as a point of a 6Z-dimensional 

statistical manifold with coordinates given by the numerical values of the expectations (x^\ =(*) 9^ and Axo = 



-(«) „ /-(«) 



constraint equations, 



= ( 2 ) 9^ . The available information can be written in the form of the following 61 information 
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The probability distributions Pi in (fT]) are constrained by the conditions of normalization, 
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Information theory identifies the Gaussian distribution as the maximum entropy distribution if only the expectation 
value and the variance are known [2l|. ME methods allow us to associate a probability distribution P (^X\Q^j to 

each point in the space of states [7j. The distribution that best reflects the information contained in the prior 
distribution m (x^j updated by the information (/xa*^ , Ax^^j is obtained by maximizing the relative entropy 
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As a working hypothesis, the prior m (x) is se t to be uniform since we assume the lack of prior available information 

about the system (postulate of equal a priori probabilities). Upon maximizing ([3]), given the constraints ([T]) and l[2|). 
we obtain 
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and, in the standard notation for Gaussians, ^8^ = f ^xl"'^) = ^8^ = f Azi"' = a^. The probability 

distribution ([4]) encodes the available information concerning the system. Note we assumed uncoupled constraints 

(a) 

among microvariables Xa . In other words, we assumed that information about correlations between the microvariables 
need not to be tracked. This assumption leads to the simplified product rule (HJ). However, coupled constraints would 
lead to a generalized product rule in (j4]) and to an information metric tensor with non-trivial off-diagonal elements 



(covariance terms). For instance, the total probability distribution P 
distributed microvariables x and y reads 



x, y\fj, x , <7 X , n v , <7y) of two dependent Gaussian 



P(X, V\Vx, a x, V>y, Vy) = 
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where r S (—1, + 1) is the correlation coefficient given by 

_ (( x ~ ( x ))(y - (v))) = i x v) - i x ) (y) 

y/{x- {x))y/{y- (y)) cr x a y 
The information metric tensor induced by © is 16], 
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where i, j = 1, 2, 3, 4. The Ricci curvature scalar associated with manifold characterized by 



is given by 



R 



8 (r 2 - 2) + 2r 2 (3r 2 - 2) 
8(r 2 -l) 



(9) 



It is clear that in the limit r — > 0, the off-diagonal elements of vanish and the scalar 1Z reduces to the result obtained 
in [l(|, namely 1Z — — 2 < 0. We could have in principle considered a correlated Gaussian process characterized by 
two correlated Gaussian microvariables x and y. Such a process would lead to an ED on a five-dimensional statistical 
manifold whose elements are probability distributions of the form P {x 1 y\[i x , a x , fi y , a y , a xy ) coordinatized by the 
expectation values y, x and \i y as well as the square-root of the three independent elements of the symmetric covariance 
matrix, namely, a x , a y and a xy . However, in view of the computational difficulty in obtaining analytical expressions 
for the elements of a 5 x 5 information metric with a xy playing the role of a macro-dynamical variable, we have chosen 
to consider the ED on a four-dimensional statistical manifold whose elements are given in (|6]) . Correlation terms may 
be fictitious. They may arise for instance from coordinate transformations. On the other hand, correlations may 
arise from external fields in which the system is immersed. In such situations, correlations among Xa effectively 
describe interaction between the microvariables and the external fields. Such generalizations would require more 
delicate analysis. 

We cannot determine the evolution of microstates of the system since the available information is insufficient. Not 
only is the information available insufficient but we also do not know the equation of motion. In fact there is no 
standard "equation of motion". Instead we can ask: how close are the two total distributions with parameters 



a[ a) ) and 



(a) 



cri"- 1 + daa*')! Once the states of the system have been defined, the next step concerns 
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the problem of quantifying the notion of change from the macrostate 8 to the macrostate O + dO. A convenient 
measure of change is distance. The measure we seek is given by the dimcnsionless distance ds between P (X |e) and 

p(x e + do), 



ds 2 = g^dO^dO" with /i, v = 1, 2,.., 61 



where 
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is the Fisher-Rao information metric. Substituting ([4]) into (fTTj) . the metric <? M „ on A4 S becomes a 6/ x 61 matrix M 
made up of 31 blocks M 2X 2 with dimension 2x2 given by, 



M 2x2 = 



,(<*) 





2 x ( ai a) 



(12) 



with a = 1, 2,...., / and a = 1,2,3. From ([TT]), the "length" element (TO]) reads, 
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(13) 



We bring attention to the fact that the metric structure of M. s is an emergent (not fundamental) structure. It arises 
only after assigning a probability distribution P 0^ to each state 9. 



III. INFORMATION-GEOMETRIC INDICATORS OF CHAOS WITHIN THE IGAC 



The relevant indicators of chaoticity within the IGAC are the Ricci scalar curvature IZm, ( or ; more correctly, the 
sectional curvature /C» s ), the Jacobi vector field intensity Jms an d the IGE <Sm s once the line element on the curved 
statistical manifold M. s underlying the entropic dynamics has been specified. 



A. Ricci Scalar Curvature, Anisotropy and Compactness 



Given the Fisher-Rao information metric, we use standard differential geometry methods applied to the space of 
probability distributions to characterize the geometric properties of M s . Recall that the Ricci scalar curvature TZ is 
given by 

TZ = g^R^, (14) 
where g^ u g vp — 5% so that g^ v = (s^ii/) - ■ The Ricci tensor is given by, 

Rtiy — djT^ — <9„r^ A + r^„r^ — rjj 7 r^. (15) 

The Christoffel symbols appearing in the Ricci tensor are defined in the standard manner as, 

T 1u = \g pa (d^gau + d^g^ - d a g^) . (16) 
Using (|12p and the definitions given above, we can show that the Ricci scalar curvature becomes 

RM S =R a a=J2 K ^ ^) = < °- 

The scalar curvature is the sum of all sectional curvatures of planes spanned by pairs of orthonormal basis elements 
{e p = 9e (p)} of the tangent space T p M s with p 6 A4 S , 

K(a,b) = - / -,a = Y(a,hr)e p , (18) 

where (e p , h?) = S a . Notice that the sectional curvatures completely determine the curvature tensor. From (JT7J) 
we conclude that A4 S is a 6^-dimensional statistical manifold of constant negative Ricci scalar curvature. A detailed 
analysis on the calculation of Christoffel connection coefficients using the ED formalism for a four-dimensional manifold 
of Gaussians can be found in flfjl. 
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It can be shown that A4 S is not a pseudosphere (maximally symmetric manifold). The first way this can be 
understood is from the fact that the Weyl Projective curvature tensor [12] (or the anisotropy tensor) Wa Vpa defined 
by 

Wfj, V pa = Rpvpa ~ ^ ^ (,9vo9pp 9vp9po) t (19) 

with n — 61 in the present case, is non- vanishing. In (|19[) . the quantity Rp Upa is the Riemann curvature tensor defined 
in the usual manner by 

R a 0p<j — d<jT a p p — d p T a p a + r a A CT r A pp — v a A P r A p a . (20) 

Considerations regarding the negativity of the Ricci curvature as a strong criterion of dynamical instability and the 
necessity of compactness of M. s in "true" chaotic dynamical systems would require additional investigation. 

The issue of symmetry of M s can alternatively be understood from consideration of the sectional curvature. In 
view of (fT8|) . the negativity of the Ricci scalar implies the existence of expanding directions in the configuration space 
manifold Ai s - Indeed, from (|17|) one may conclude that negative principal curvatures (extrema of sectional curvatures) 
dominate over positive ones. Thus, the negativity of the Ricci scalar is only a sufficient (not necessary) condition 
for local instability of geodesic flow. For this reason, the negativity of the scalar provides a strong criterion of local 
instability. Scenarios may arise where negative sectional curvatures are present, but the positive ones could prevail in 
the sum so that the Ricci scalar is non-negative despite the instability in the flow in those directions. Consequently, 
the signs of the sectional curvatures are of primary significance for the proper characterization of chaos. 

Yet another useful way to understand the anisotropy of the M s is the following. It is known that in n dimensions, 
there are at most independent Killing vectors (directions of symmetry of the manifold). Since M s is not a 

pseudosphere, the information metric tensor docs not admit the maximum number of Killing vectors K v defined as 



where Dp, defined as 



Ck9p„ = DpK v + D v Kp = 0, (21) 



DpK v = dpK v - TP^Kp (22) 



is the covariant derivative operator with respect to the connection T defined in (I16|) . The Lie derivative Hk9pv of the 
tensor field gp V along a given direction K measures the intrinsic variation of the field along that direction (that is, 
the metric tensor is Lie transported along the Killing vector) [23[ . Locally, a maximally symmetric space of Euclidean 
signature is either a plane, a sphere, or a hyperboloid, depending on the sign of R. In our case, none of these scenarios 
occur. As will be seen in what follows, this fact has a significant impact on the integration of the geodesic deviation 
equation on M s . At this juncture, we emphasize it is known that the anisotropy of the manifold underlying system 
dynamics plays a crucial role in the mechanism of instability. In particular, fluctuating sectional curvatures require 
also that the manifold be anisotropic. However, the connection between curvature variations along geodesies and 
anisotropy is far from clear and is currently under investigation. 

Krylov was the first to emphasize |24[ the use of 7Z < as an instability criterion in the context of an iV-body 
system (a gas) interacting via Van der Waals forces, with the ultimate hope to understand the relaxation process in a 
gas. However, Krylov neglected the problem of compactness of the configuration space manifold which is important 
for making inferences about exponential mixing of geodesic flows [251 ]. Why is compactness so significant in the 
characterization of chaos? True chaos should be identified by the occurrence of two crucial features: 1) strong 
dependence on initial conditions and exponential divergence of the Jacobi vector field intensity, i.e., stretching of 
dynamical trajectories; 2) compactness of the configuration space manifold, i.e., folding of dynamical trajectories. 
Compactness [U, is required in order to discard trivial exponential growths due to the unboundedness of the 
"volume" available to the dynamical system. In other words, the folding is necessary to have a dynamics actually 
able to mix the trajectories, making practically impossible, after a finite interval of time, to discriminate between 
trajectories which were very nearby each other at the initial time. When the space is not compact, even in presence 
of strong dependence on initial conditions, it could be possible in some instances (though not always), to distinguish 
among different trajectories originating within a small distance and then evolved subject to exponential instability. 

As a final remark, we emphasize that it is known from IG Q that there is a one-to-one relation between elements 
of the statistical manifold and the parameter space. More precisely, the statistical manifold M s is homeomorphic to 
the parameter space T>q. This implies the existence of a continuous, bijective map Hm 3 , p e , 



v e :M S 3P [X 6 U96P e (23) 
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where hjj^ Vq ^O^ — P \X 0J. The inverse image hjj^ Vq is the so-called homeomorphism map. In addi- 
tion, since homeomorphisms preserve compactness, it is sufficient to restrict ourselves to a compact subspace of the 
parameter space Vq in order to ensure that Ms is itself compact. 

B. Canonical Formalism 

The geometrization of a Hamiltonian system by transforming it to a geodesic flow is a well-known technique of 
classical mechanics associated with the name of Jacobi [1 81 ] . Transformation to geodesic motion is obtained in two 
steps: 1) conformal transformation of the metric; 2) rescaling of the time parameter [27[. The reformulation of 
dynamics in terms of a geodesic problem allows the application of a wide range of well-known geometrical techniques 
in the investigation of the solution space and properties of equations of motions. The power of the Jacobi reformulation 
is that all of the dynamical information is collected into a single geometric object - the manifold on which geodesic 
flow is induced - in which all the available manifest symmetries are retained. For instance, integrability of the system 
is connected with the existence of Killing vectors and tensors on this manifold [28|, H^] . 

In this Section we study the trajectories of the system on Ai s . We emphasize ED can be derived from a standard 
principle of least action (of Maupertuis-Euler-Lagrange- Jacobi type) [H, [3(|. The main differences are that the 
dynamics being considered here, namely ED, is defined on a space of probability distributions Ai s , not on an ordinary 
linear space V and the standard coordinates of the system are replaced by statistical macrovariables 8 P . The 
geodesic equations for the macrovariables of the Gaussian ED model are given by, 

d 2 Q» „„ dO» dQP n 

with fi — 1, 2,..., 6/. Observe the geodesic equations are nonlinear second order coupled ordinary differential equations. 
They describe a reversible dynamics whose solution is the trajectory between an initial and a final macrostate. The 
trajectory can be equally well traversed in both directions. 

1, Geodesies on M s 

We determine the explicit form of (|24p for the pairs of statistical coordinates (p,a"\ o"o )• Substituting the expres- 
sion of the Christoffel connection coefficients into (|24[) . the geodesic equations for the macrovariables and 
associated to the microstate Xa become, 

d 2 ^ 2 dfj,^ da™ _ rfVj a) 1 (da^V 1 (d^\ 2 _ 
dr* a W dr dr U ' dr> ff (°0 { dr ) + 2c7 M [dr ) [ 

with a = 1, 2,...., I and a = 1, 2, 3. This is a set of coupled ordinary differential equations, whose solutions are 



K°>) 

M w (r) = H£T 



-sinh 



2" I 

(26) 



(a) , x _ Bj°> cxp(-/3i,°)r) ( Q ) 

y T > (W>y~ ' 



The quantities B^"\ ci°\ f3^ are real integration constants and they can be evaluated once the boundary conditions 
are specified. We observe that since every geodesic is well-defined for all temporal parameters t, M. s constitutes a 
geodesically complete manifold [3lj . It is therefore a natural setting within which one may consider global questions 

< +oo Vt e R+, 



and search for a weak criterion of chaos |2|. Furthermore, since /Lti"" 1 (t) < +oo and ct^' (t) 

Va = 1, 2, 3 and Va = 1,.., JV", the parameter space |©| (homeomorphic to M s ) is compact. The compactness of the 

configuration space manifold Ai s assures the folding mechanism of information-dynamical trajectories (the folding 
mechanism is a key-feature of true chaos, 0). 
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It is known (30j that the Riemannian curvature of a manifold is intimately related to the behavior of geodesies on it. 
If the Riemannian curvature of a manifold is negative, geodesies (initially parallel) rapidly diverge from one another. 

For the sake of simplicity, we assume very special initial conditions: Ba = S, (3^ = A e M + , Ca = 0, Vet = 1, 
2,...., I and a = 1, 2, 3. However, the conclusions drawn can be generalized to more arbitrary initial conditions. We 
observe that since every maximal geodesic is well-defined for all temporal parameters r, M. s constitute a geodesically 
complete manifold [3l|. It is therefore a natural setting within which one may consider global questions and search 
for a weak criterion of chaos [2] • 



C. Exponential divergence of the Jacobi field intensity 



The actual interest of the Riemannian formulation of the dynamics stems form the possibility of studying the 
instability of natural motions through the instability of geodesies of a suitable manifold, a circumstance that has 
several advantages. First of all a powerful mathematical tool exists to investigate the stability or instability of a 
geodesic flow: the Jacobi-Levi-Civita equation for geodesic spread [32j]. The JLC-equation describes covariantly how 
nearby geodesies locally scatter. It is a familiar object both in Riemannian geometry and theoretical physics (it 
is of fundamental interest in experimental General Relativity). Moreover the JLC-equation relates the stability or 
instability of a geodesic flow with curvature properties of the ambient manifold, thus opening a wide and largely 
unexplored field of investigation of the connections among geometry, topology and geodesic instability, hence chaos. 

Consider the behavior of the one-parameter family of neighboring geodesies Tg m , (A) = {©xV ( T > 
where 



AGR+ 



^ (r; A) 



—2 



1 



2A cosh (2At) - sinh (2At) 



°i a) (r; A) = H 



cosh (2At) - sinh (2Ar) 



SA 2 



(27) 



cosh (At) — sinh (At) 



,xA- 



with a = 1, 2,...., I and a = 1, 2, 3. The relative geodesic spread on a (non-maximally symmetric) curved manifold as 
M s is characterized by the Jacobi-Levi-Civita equation, the natural tool to tackle dynamical chaos [23l. |32T|, 



(28) 



where the covariant derivative 



D 2 S6» 



Dt 2 



d 2 5Q» 
dr 2 

n R L r 



is defined as [3 

dse a d® 13 



. ore 

^ dr dr 

dQ a d® 13 



af} pa dr dr 
and the Jacobi vector field is given by [34| . 



«? dU dr 2 + a ^ dr dr 



( d<3» (t; A) 



SX. 



(29) 



(30) 



Notice that the JLC-equation appears intractable already at rather small I. For isotropic manifolds, the JLC-equation 
can be reduced to the simple form [32l |. 



D 2 J^ 
Dt 2 



KJ^ =0,(i = 1,. 



61 



(31) 



where K is the constant value assumed throughout the manifold by the sectional curvature. The sectional curvature 
of manifold A4 S is the 6Z-dimensional generalization of the Gaussian curvature of two-dimensional surfaces of K 3 . If 
K < 0, unstable solutions of the equation (|31[) are of the form 



J(t) 



=uj (0) sinh i 



(32) 
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once the initial conditions are assigned as J (0) = 0, dJ }^ = uj (0) and K < 0. Equation (f28|) forms a system of 61 
coupled ordinary differential equations linear in the components of the deviation vector held (|30|) but nonlinear in 
derivatives of the metric (fTTj) . It describes the linearized geodesic flow: the linearization ignores the relative velocity of 
the geodesies. When the geodesies are neighboring but their relative velocit y is arbitrary, the corresponding nonlinear 
geodesic deviation equation is the so-called generalized Jacobi equation |35l . l3a |. The nonlinearity is due to the 
existence of velocity-dependent terms in the system. Neighboring geodesies accelerate relative to each other with a 
rate directly measured by the curvature tensor R a p~/6- Substituting (|2T[) in (|28[) and neglecting the exponentially 
decaying terms in SQ^ and its derivatives, integration of (f2"5)) leads to the following asymptotic exponential growth 
of the Jacobi vector field intensity (a classical feature of chaos), 



Jm s = \\J\\ = {g^JTY 

Finally, we point out that in our approach the quantity Ay, 



T 3le^ 



(33) 



del' 



f 



A / = lim — In 





(r) 




(0) 



(34) 



would play the role of the conventional Lyapunov exponents. 



IV. LINEARITY OF THE INFORMATION GEOMETRODYNAMICAL ENTROPY 

We investigate the stability of the trajectories of the ED model considered on M s . It is known [3(| that the 
Riemannian curvature of a manifold is closely connected with the behavior of the geodesies on it. If the Riemannian 
curvature of a manifold is negative, geodesies (initially parallel) rapidly diverge from one another. For the sake of 
simplicity, we assume very special initial conditions: = S, (3^ = A e R+, = 0, V a = 1, 2,...., / and 

a = 1, 2, 3 . However, the conclusion we reach can be generalized to more arbitrary initial conditions. Recall M. s is 

the space of probability distributions |p 9^ j labeled by 61 statistical parameters O. These parameters arc the 
coordinates for the point P, and in these coordinates a volume element cIVm 3 reads, 

dV Ms = V?d 6 'e = n fl -^Tld^da^. (35) 
o=i«=i 

The volume of an extended region AFa/U ( r ! ^) °f M s is defined by, 

A^(r;A)^niI /, , L J^d^da^ (36) 



where Ha and (r) are given in (|27p and where the scalar A is the chosen quantity used to define the one- 
parameter family of geodesies Tg m . W ^ {®!m s ( r > ^)}AeR+' 6 '' ( l uan *i*y tnat encodes relevant information 
about the stability of neighboring volume elements is the the average volume Vm 3 (t; A), 

r 

V Ms (r; A) = (AV M . (r'; A)) r - / AV Ms (r'; A) dr' ™ . (37) 



T 







The ratio with Vm 3 (t) in ([37| representing the temporal average of the 3/-fold integral over trajectories of 

maximum probability (geodesies) is a measure of the number of the accessible macrostates in configuration (statistical) 
manifold M. s after a finite temporal increment r. In other words, Vm s (t) can be interpreted as the temporal evolution 
of the system's uncertainty volume Vm s (0). For instance Vm s (0) ma y be a spherical volume of initial points whose 
center is a given point on the attractor and whose surface consists of configuration points from nearby trajectories. An 
attractor is a subset of the manifold M. s toward which almost all sufficiently close trajectories converge asymptotically, 
covering it densely as the time goes on. Strange attractors are called chaotic attractors. Chaotic attractors have at 
least one finite positive Lyapunov exponent [37j |. As the center of Vm s (0) an( i its surface points evolve in time, the 
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spherical volume becomes an ellipsoid with principal axes in the directions of contraction and expansion. The average 
rates of expansion and contraction along the principal axes are the Lyapunov exponents [38j . 

The asymptotic regime of diffusive evolution in (|37|) describes the exponential increase of average volume elements 
onM s . The exponential instability characteristic of chaos forces the system to rapidly explore large areas (volumes) 
of the statistical manifolds. From equation (|37p . we notice that the parameter A characterizes the exponential growth 
rate of average statistical volumes Vm 3 (t; A) in M s . This suggests that A may play the same role ordinarily played 
by Lyapunov exponents [3^ |. It is interesting to note that this asymptotic behavior appears also in the conventional 
description of quantum chaos where the von Neumann entropy increases linearly at a rate determined by the Lyapunov 
exponents. The linear increase of entropy as a quantum chaos criterion was introduced by Zurek and Paz [4|. In our 
information-geometric approach a relevant quantity that can be useful to study the degree of instability characterizing 
the ED model is the information-geometrodynamical entropy (IGE) defined as [l(| , 



S 



dcf 



Ms 



lim logV.M s (r; A) 



(38) 



The IGE is intended to capture the temporal complexity (chaoticity) of ED theoretical models on curved statistical 
manifolds A4 S by considering the asymptotic temporal behaviors of the average statistical volumes occupied by the 
evolving macrovariables labelling points on M s - Substituting (|37|) in (|38|) . we obtain 



?M, = km log 

r — >oo | T 



I 3 

nn 

a=l o=l J ^ 



>(0) 



dr' 



31Xt. 



(39) 



Before discussing the meaning of ([39]) . recall that in conventional approaches to chaos the notion of entropy is 
introduced, in both classical and quantum physics, as the missing information about the systems fine-grained state 
@, El- For a classical system, suppose that the phase space is partitioned into very fine-grained cells of uniform 
volume Aw, labelled by an index j. If one does not know which cell the system occupies, one assigns probabilities pj 
to the various cells; equivalently, in the limit of infinitesimal cells, one can use a phase-space density p(Xj) — J^j. 
Then, in a classical chaotic evolution, the asymptotic expression of the information needed to characterize a particular 
coarse-grained trajectory out to time r is given by the Shannon information entropy (measured in bits), 



S. 



(chaotic) 
classical 



dXp (X)log 2 ( P (X)Av) 



(40) 



where p (X) is the phase-space density and pj — is the probability for the corresponding coarse-grained trajectory. 

is the missing information about which fine-grained cell the system occupies. The quantity IC represents 



S, 



(chaotic) 
classical 



the linear rate of information increase and it is called the Kolmogorov-Sinai entropy (or metric entropy) {IC is the 
sum of positive Lyapunov exponents, IC — Y]j Xj ). IC quantifies the degree of classical chaos. The Kolmogorov- 
Sinai entropy provides a measure of the rate at which information is lost by an evolving chaotic system ( IC has 
dimension entropy/time) and has its roots in the definition of the Shannon entropy. It is worthwhile emphasizing 
that the quantity that grows asymptotically as Kt is really the average of the information on the left side of equation 
(|40p . This distinction can however be ignored provided we assume the chaotic system has roughly constant Lyapunov 
exponents over the accessible region of phase space. In quantum mechanics the fine-grained alternatives are normalized 
state vectors in Hilbert space. From a set of probabilities for various state vectors, one can construct a density operator 



(41) 



The normalization of the density operator, tr (p) = 1, implies that the eigenvalues make up a normalized probability 
distribution. The von Neumann entro py (natural generalization of both Boltzmann's and Shannon's entropy) of the 
density operator p (measured in bits) [4l| , 



-.(chaotic) 
^quantum 



= -tr (plog 2 p) = - V Xj log 2 A 



ICqT 



(42) 



can be thought of as the missing information about which eigenvector the system is in. Entropy quantifies the degree 
of unpredictability about the system's fine-grained state. In quantum mechanics, the von Neumann entropy plays a 
role analogous to that played by the Shannon entropy in classical probability theory. They are both functionals of 
the state, are both monotone under a relevant kind of mapping, and can both be singled out uniquely by natural 
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requirements, von Neumann's entropy reduces to the Shannon entropy for diagonal density matrices. However, in 
general the von Neumann entropy is a subtler object than its classical counterpart. The quantity tC q in (|42[) can be 
interpreted as the non-commutative (quantum theory is a non-commutative probability theory) quantum analog of the 
Kolmogorov-Sinai dynamical entropy, the so-called quantum dynamical entropy [43[ . Examples of quantum dynamical 
entropies applied to quantum chaos and quantum information theory are the Alicki-Fannes (AF) [44j entropy and 
the Connes-Narnhofer-Thirring (CNT) [45[ entropy. Both the AF and CNT entropy coincide with the KS entropy on 
classical dynamical systems. They also coincide on finite-dimensional quantum systems. However, they differ when 
moving from finite to infinite quantum systems. 

Recall that decoherence is the loss of phase coherence between the set of preferred quantum states in the Hilbert 
space of the system due to the interaction with the environment. Moreover, decoherence induces transitions from 
quantum to classical systems. Therefore, classicality is an emergent property of an open quantum system. Motivated 
by such considerations, Zurek and Paz investigated implications of the process of decoherence for quantum chaos. 
They considered a chaotic system, a single unstable harmonic oscillator characterized by a potential V (x) = — ^f- 
(A is the Lyapunov exponent), coupled to an external environment. In the reversible classical limit [42j, the von 
Neumann entropy of such a system increases linearly at a rate determined by the Lyapunov exponent, 

45™ (Zurek-Paz) T ~°° At. (43) 

Notice that the consideration of 3/ uncoupled identical unstable harmonic oscillators characterized by potentials 
Vi (x) = — ^f— (\i = Aj-; i, j = 1, 2,..., 3/) would simply lead to 

S<S£ (Zurek-Paz) MXr. (44) 

The resemblance of equations (|3T)1) and (|4"4"]) is remarkable and a more detailed discussion about this analogy is 
presented in (l6j where an information-geometric analogue of the Zurek-Paz quantum chaos criterion in the classical 
reversible limit is proposed [J] . This analogy is illustrated applying the IGAC to a set of n-uncoupled three-dimensional 
anisotropic inverted harmonic oscillators characterized by a Ohmic distributed frequency spectrum. 

The entropy- like quantity Sm , m (|39P is the asymptotic limit of the natural logarithm of the statistical weight 
(AVm 3 ) t defined on Ai s and it grows linearly in time, a quantum feature of chaos. Indeed, equation (|39[) may be 
considered the information-geometric analog of the Zurek-Paz chaos criterion. In our chaotic ED Gaussian model, the 
IGE production rate is determined by the information-geometric parameter A characterizing the exponential growth 
rate of average statistical volumes Vm s ( t ! ^) m M s - 

In conclusion, for the example under investigation, we have 

T — *00 T — > OO \ 

H Ms = -3i, S Ms « 3ZAt, J Ms ps 3Ze Ar . (45) 

The IGE grows linearly as a function of the number of Gaussian-distributed microstates of the system. This supports 
the fact that Sm, may be a useful measure of temporal complexity [461 ]. Furthermore, these three indicators of chaotic- 
ity, the Ricci scalar curvature TZm s j the information-geometric entropy Sm 3 and the Jacobi vector field intensity Jm s 
are proportional to 3^, the dimension of the microspace with microstates |a| underlying our chaotic ED Gaussian 

model. This proportionality leads to the conclusion that there is a substantial link among these information-geometric 
measures of chaoticity since they are all extensive functions of the dimensionality of the microspace underlying the 
macroscopic chaotic entropic dynamics (see (|45[1 ). Curvature, information-geometrodynamical entropy and Jacobi 
field intensity are linked within our formalism. We are aware that our findings are reliable in the restrictive assump- 
tion of Gaussianity. However, we believe that with some additional technical machinery, more general conclusions can 
be achieved and this connection among indicators of chaoticity may be strengthened. 



V. INFORMATION GEOMETRY OF QUANTUM ENERGY LEVEL STATISTICS: AN APPLICATION 

TO ENTANGLEMENT IN QUANTUM SPIN CHAINS 

In what follows, we apply the IGAC to study the entropic dynamics on curved statistical manifolds induced by 
classical probability distributions of common use in the study of regular and chaotic quantum energy level statistics. 
In doing so, we suggest an information-geometric characterization of a special class of regular and chaotic quantum 
energy level statistics. More precisely, we present an information-geometric analogue of the logarithmic and linear 
entanglement entropy growth in regular and quantum chaotic spin chains, respectively. 
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A. The Information Geometry of the Poisson and Wigner-Dyson Distributions 

The theory of quantum chaos (quantum mechanics of systems whose classical dynamics are chaotic) is not primarily 
related to few-body physics. Indeed, in real physical systems such as many-electron atoms and heavy nuclei, the 
origin of complex behavior is the quite strong interaction among many particles. To deal with such systems, a famous 
statistical approach has been developed which is based upon the Random Matrix Theory (RMT) j4?|]. The main 
idea of this approach is to neglect the detailed description of the motion and to treat these systems statistically 
bearing in mind that the interaction among particles is so complex and strong that generic properties are expected 
to emerge. The simplest models of RMT are full random matrices of a given symmetry. One of the main results of 
RMT is the prediction of a specific kind of correlations of the energy spectra of complex quantum systems. Among 
many characteristics of these correlations, the most popular one is the distribution of spacings between nearest energy 
levels in the spectra. The exact analytical expression of this distribution is very complicated; instead, one uses 
the so-called Wigner-Dyson surmise (a very simple expression which gives a very good approximation to the exact 
result). The known manifestation of quantum chaos is the so-called Wigner-Dyson (WD) distribution for spacings 
between neighboring levels in the spectrum. In the other limiting case of completely integrable (regular) systems, the 
distribution turns out to be very close to the Poissonian one. A distinctive property of the WD distribution is the 
repulsion between neighboring levels in the spectra; the degree of this repulsion (linear, quadratic or quartic) depends 
on the symmetry of random matrices. For systems without time reversal invariance the relevant ensemble of random 
matrices is the Gaussian Unitary Ensemble (GUE) [53], characterized by the probability distribution 

32 / 4 \ 

Pgue (0) = — 2 exp I 9 2 ) , (quadratic repulsion) (46) 

n z \ 7T ) 

where 9 is the average spacing of the energy levels. For systems invariant with respect to time reversal the ensemble 
is the Gaussian Orthogonal Ensemble (GOE) [47] . 

Pgoe (0) — ^-6>exp (~^ 2 J > (linear repulsion). (47) 

For systems with time reversal invariance but with half-integer spin, the energy is described by the Gaussian Symplectic 
Ensemble (GSE) of random matrices (47| . 

" JS n4 ( 64, 



Pgse (0) = t^^O ex P ( J ' (q uartic repulsion). (48) 

Equations (|4"6"| . (|4T)l and (|4"8"| are standard accepted conjectures. Besides energy level statistics in the extreme 
integrable (Poisson) and chaotic (Wigner-Dyson) regimes, there is also energy level statistics in the mixed regime, i.e., 
such having a mixed classical dynamics where regular and chaotic regions coexist in the phase space. A convenient 
and often successful paramctrization of the correct probability distribution in the transition region between Poisson 
and WD distributions is provided by the Brody interpolation formula [48| , 

pjBrody) W = 7 ^+l)exp(- 7 ^ +1 ). (49) 

r r/3+2l "I 13+1 

where 7 = < T > and T ((3) is the Euler Gamma function. This distribution is normalized and, by construc- 

tion, has mean spacing (6) = 1. We recover the Poisson case by taking (3 = while the Wigner case is recovered for 
(3=1. However, a criticism of the Brody distribution is the lack of a first principles justification for its validity. The 
fact remains that it does fit the specific results found when considering explicit model systems. It is essentially an 
ad hoc one-parameter family of distributions and has no deep physical background, but it does interpolate between 
Poisson and Wigner-Dyson in a simple, effective manner. Our objective here is to apply our information-geometric 
formalism (based on statistical inference methods) to Wigner-Dyson and Poisson probability distributions. 

Most of the probability distributions arise from the maximum entropy formalism as a result of some simple state- 
ments concerning averages. Not all distribution are generated in this way. Some distributions are generated by 
combining the results of simple cases (multinomial from a binomial) . Other distributions are found as a result of a 
change of variable (Cauchy distribution). For instance, the Weibull distribution 49] can be obtained from an expo- 
nential distribution as a result of a power law transformation. Assume our knowledge of the microstate x is encoded 
in an exponential distribution, 

p{x\ff) = \e-*, (50) 
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where x may be considered the spacing of the energy levels while 9 is the average spacing, 6 — (x) . Note that the study 
of probability distributions could, in principle, be restricted to the exponential type since an arbitrary distribution 
can be represented in exponential form. It is said that the exponential family of distributions is dense in the totality of 
probability distributions (50j . We can re-express x £ X in p (x\9) in terms of another random variable y = f (x) G y, 
assuming / is an invertible mapping. For instance, consider the power law transformation 



x^y = f(x) 



(51) 



We clearly have, 



x 



Vo\d{x) -> p nC w(y)= / dxpou (x) 6 (y - f (x)) 



S (r 1 (y) - x) 



dxp \d (x) 



X 



Po\d (x) 



(52) 



Therefore, considering and ([5T|) . equation ([52")) leads to 



It is worthwhile emphasizing that since 



P„cw(y) =n-e ° v "y n x . 



does not depend on 9 and since Jdy = jdx 

y x 



(53) 



, we have 



/ 



dyf>ncw (y) <9 M log 

Pncw (y) d v log (y) = / dyp old (x) <9 M logp id {x) 8 V \ogp Q \ A (x) 



(54) 



x 



Equation (f5"4")) leads to conclude that the Fisher-Rao information metric g^ v is invariant under transformations of 
the random variable. For the sake of completeness, let us show that the information metric is also covariant under 

reparametrization. Suppose that f^ij * s a new se * °f coordinates, specified in terms of the old set through the 
invertible relationship 9^ — 9^ (9). Defining p. (x) = P s( *s ( x )i we are then able to compute the new metric tensor 
v (e\ in terms of (9). Indeed, since -^fP 6 = §p m*P ($y we obtain 

>6> p 89° 



97 



-(») 



89 89 



8=6(8) 



(55) 



Letting g = jrr, from (|53p we obtain the Weibull probability distribution, 



A 



Pwcibuii (y|A) = j (j^ 

Moreover, letting n = 2, y — A and A = from ([56]) we obtain the standard Wigner-Dyson distribution, 

7rA 



PWigncr- Dyson 

(A\D) 



2D 2 ' 2 



(56) 



(57) 



In conventional notations, A is the spacing between two neighboring energy levels and D is the average spacing [5l| . 
Recall that the Fisher-Rao information metric G^J (9) of a Poissonian probability distribution p(x\9) is defined as, 



G$ (9) = I dxp (x\9) 8^ \ogp (x\9) d v \ogp (x\9) with d„ = ^ 



where p (x\9) is given by, 



p(xh 



exp 



(58) 



(59) 
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The Poisson line element (ds 2 )^ . is defined as, 

\ / Poisson 7 

(ds 2 ) p . = G$ (e) dO^dO" = \d8 2 . (60) 

V /Poisson ftv \ ) qI V / 

The Fisher-Rao information metric G^v ((/>) of a Wigner-Dyson probability distribution q (y\(j)) is defined as, 

C Pi 

G^ D) {4>) = J dyqiy^d^ogqiy^dAomiv^) with = — , (61) 



where q(y\4>) is given by, 



Notice that q(y\<j>) in (j6"2|) is equivalent to pwigncr-Dyson (A\D) in ([57]) with y = A and <j> = D. The Wigner-Dyson 
line element ((is 2 ),,-. _ is defined as, 

V / Wigner-Dyson ' 

(^ 2 )wi g ne r - DyS o„=^ D) W^^- (63) 

Notice that the Poisson distribution and the Wigner-Dyson distributions are related through the combination of a 
change of random variable and a new reparametrization, namely 

1 {y\4>) =p{x{y)\8 ((f))) J (y) (64) 

where, 



(y) = \y 2 ,6(<f>) = ^\, J(y) = 



dx (y) 



(65) 



dy 

Considering equations (f54|) and (|55| . the Wigner-Dyson line element (^ s2 ) wigncr Dyson b ecomes 

(^ 2 )w Ig „e,Dyso„ = WW 

= G<$ [9 m dd dd [4>T = ^dtf. (66) 

Equations (|60|) and (|66l) will be used in our IGAC in quantum spin chains systems. Before considering such 
information-geometric characterization of quantum energy level statistics for regular and chaotic spin chains im- 
mersed in an external magnetic field, we briefly review the main points of the more standard approach to these 
topics. 



B. Entanglement in quantum spin chains: standard formalism 



One of the most important concepts in quantum information theory is that of entanglement, an intrinsic property 
of composite quantum systems. Entanglement plays an essential role in many-body quantum phenomena, such as 
superconductivity [52j and quantum phase transitions [53j ] . Moreover, it is an important concept in quantum compu- 
tation and information processing [54j |. An excellent theoretical framework for investigating entanglement properties 
is offered by spin chains. Quantum spin chains belong to the most studied models of quantum statistical mechanics. 
However, only for a few types of models have the thermal and ground state structures have been determined. This 
is mainly a consequence of the complicated correlations that can arise among quantum states. These strong corre- 
lations can even be present in pure quantum states, while classical pure states can only have a trivial product state 
structure. Unlike the classical case, the restrictions of pure states on the quantum spin chain to local subsystems are 
typically mixed states. This type of correlation between subsystems is commonly referred to as entanglement. The 
von Neumann entropy, defined as, 

Svon Neumann = ~tr (plogp) , (67) 

is a standard measure of the nonpurity of the reduced density matrix p, thus it is a very useful quantity in the 
description of entanglement (55j . Several simple models of spin chains can be studied analytically and there also exist 
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efficient numerical techniques. There are two widely used methods of characterizing entanglement in spin chains. 
The first of these describes the entanglement between two spins in the chain with a quantity called concurrence [5(| . 
The other one measures entanglement of a block of spins with the rest of the chain with the von Neumann entropy 
when the chain is in its ground state [57j ] . As a side remark, we emphasize that entanglement entropy does not refer 
exclusively to the characterization of quantum systems in their ground states, but, more generally, it refers to any 
many particles quantum dynamical states undergoing a unitary time evolution. The method used in (57| is known as 
the density matrix renormalization group method (DMRG) [58| . It is based on the fact that many degrees of freedom 
are redundant in quantum state description; therefore, the system is adequately described by taking into account 
maximally entangled components only, von Neumann entropy is supposed to play an important role in quantifying 
the essential subspace of a reduced density matrix. The possibility of compressing such density matrices from its 
full dimension to a much smaller subspace without significant loss of information is the starting point of the DMRG 
analysis. Classical complexity of quantum states can be characterized by a mixed state entanglement entropy. For 
instance, the von Neumann entropy of a block of L neighboring spins in a XX chain, describing entanglement of the 
block with the rest of the chain is given by (BUI ]. 

Sl = -tr (p L \ogp L ) L «°° logL. (68) 

The reduced density matrix p L is obtained from the ground state \^ g ) of the chain by tracing out external degrees of 
freedom, 

p L = tr N _ L (v^l , H |* fl ) = E g \* B ) . (69) 
The Hamiltonian U in J69J is given bv [57115^. 

JV N 

« = -£ ( s f s t+i + s i s hi) fr£*f. ( 7 °) 
i=i i=i 

where sf (a = x, y, z) are the Pauli spin matrices at sites I — 1, 2,.., N of a periodic chain and h is the magnetic 
field. The logarithmic growth of the entanglement entropy is a general consequence of the fact that in one dimen- 
sional systems near quantum phase transitions, the entropy is a logarithmic function of the size of the system [60| . 
Furthermore, there is also a time dependent version of this DMRG, known as r-DMRG [61]. This method is used 
to study the evolution of pure states, density matrices and operators. Note that the classical complexity of quantum 
operators can be characterized using the operator space entanglement entropy of a density operator, not the state 
entanglement entropy of a mixed state. For the evolution of density matrices and operators, a superket corresponding 
to an operator O expanded in the computational basis of products of local operators is considered. For instance, for 
a chain of n-qubits, a basis of 4™ Pauli operators is used, a sa ® ... <£> a 3 " -1 , with Sj £ {0, x, y, z} and <r° = I. The 
key idea of r-DMRG is to represent any operator in a matrix product form [621 ] , 

O = J> (Aq ...A s „Si) a s ° ® .... <g> (71) 

in terms of the 4n matrices Aj 3 of fixed dimension D. The number of parameters in the matrix product (MPO) 
representation of the operator is 4nD 2 . The minimal D required, D £ (r), is equal to the maximal rank of the reduced 
super density matrix over bipartitions of the chain. The way entropy can be computed in the spaces of operators 
can be found in [62| . The r-DMRG method is very efficient in classical simulations of many body quantum dynamics 
requiring that the computational costs grow polynomially in time and, consequently, that the entanglement entropy 
grows no faster than logarithmically. However, it is known that the asymptotic behavior of computational costs and 
entanglement entropies of integrable and chaotic Ising spin chains are very different [l7T ] . Here Prosen considered the 
question of time efficiency implementing an up-to-date version of the r-DMRG for a family of Ising spin ~ chains 
in arbitrary oriented magnetic field, which undergoes a transition from integrable (transverse Ising) to nonintegrablc 
chaotic regime as the magnetic field is varied. An integrable (regular) Ising chain in a general homogeneous transverse 
magnetic field is defined through the Hamiltonian 7^( re s ular ) = 7Y (0, 2), where 

n — 2 n— 1 

h (k, hy) = + E 0» x °i + • ( 72 ) 

3=0 J=0 

In this case, the computational cost shows a polynomial growth in time, £)^ rcgulai ) ( r ) r (x°° r, while the entanglement 
entropy is characterized by logarithmic growth, 

5 (rcgular) r «°°c log r + c'. (73) 



16 



The constant c depends exclusively on the value of the fixed transverse magnetic field intensity B± , while d depends 
on B± and on the choice of the initial local operators of finite index used to calculate the operator space entanglement 
entropy. Instead, a quantum chaotic Ising chain in a general homogeneous tilted magnetic field is defined through 
the Hamiltonian H (chaotic) = H(1,1), where U is defined in ([72]). In this case, the computational cost shows an 
exponential growth in time, £)( chaotlc ) ( T ) T ^ ex p (JC q r) while the entanglement entropy is characterized by linear 
growth, 

^(chaotic) = ™ ^ ^ 

The quantity K, q is a constant, asymptotically independent of the number of indexes of the initial local operators used 
to calculate the operator space entropy, that depends only on the Hamiltonian evolution and not on the details of the 
initial state observable or error measures, and can be interpreted as a kind of quantum dynamical entropy. 

It is well known the quantum description of chaos is characterized by a radical change in the statistics of quantum 
energy levels (63|. The transition to chaos in the classical case is associated with a drastic change in the statistics 
of the nearest-neighbor spacings of quantum energy levels. In the regular regime the distribution agrees with the 
Poisson statistics while in the chaotic regime the Wigner-Dyson distribution works very well. Uncorrelated energy 
levels are characteristic of quantum systems corresponding to a classically regular motion while a level repulsion (a 
suppression of small energy level spacing) is typical for systems which are classically chaotic. A standard quantum 
example is provided by the study of energy level statistics of an Hydrogen atom in a strong magnetic field. It is known 
that level spacing distribution (LSD) is a standard indicator of quantum chaos [64 |. It displays characteristic level 
repulsion for strongly nonintegrable quantum systems, whereas for integrable systems there is no repulsion due to 
existence of conservation laws and quantum numbers. In [l7|, the authors calculate the LSD of the spectra of 7^( re s ular ) 
and H^ haotic \ They find that for H^ ceulal \ the nearest neighbor LSD is described by a Poisson distribution. For 
H (chaotic) , they find the nearest neighbor LSD is described by a Wigner-Dyson distribution. Therefore, they conclude 
that H (regular) and ft (chaotic) indeed represent generic regular and quantum chaotic systems, respectively. 

In the next paragraph, we will encode the relevant information about the spin-chain in a suitable composite- 
probability distribution taking into account the quantum spin chains and the configurations of the external magnetic 
field in which they are immersed. 



VI. AN INFORMATION GEOMETRIC MODEL OF REGULAR AND CHAOTIC QUANTUM SPIN 

CHAINS 

A. Integrable Statistical Model: Poisson coupled to an Exponential Bath 

Recall that in the ME method Q , the selection of relevant variables is made on the basis of intuition guided by 
experiment; it is essentially a matter of trial and error. The variables should include those that can be controlled 
or experimentally observed, but there are cases where others must also be considered. Our objective here is to 
choose the relevant microvariables of the system and select the relevant information concerning each one of them. In 
the integrable case, the Hamiltonian 7^( re s ular ) describes an antiferromagnetic Ising chain immersed in a transverse 
homogeneous magnetic field ^transverse = B± B ± with the level spacing distribution of its spectrum given by the 
Poisson distribution 

(Poisson) / I \ 1 / \ fnr-\ 

Pa C^aIMa) = — ex P : ( 75 ) 

Ma V Ma / 

where the microvariable xa is the spacing of the energy levels and the macrovariable /i A is the average spacing. The 
chain is immersed in the transverse magnetic field which has just one component B± in the Hamiltonian 7Y( re s ular ). 
We translate this piece of information in our IGA formalism, coupling the probability (|75j) to an exponential bath 

(exponential) / i \ • 1 

Pb C^bIMb) g lven b y 

(exponential) / i \ 1 / \ t^n\ 

Pb (^b Mb) = — exp , 76) 

Mb V Mb/ 

where the microvariable xb is the intensity of the magnetic field and the macrovariable /i B is the average intensity. 
More correctly, xb should be the energy arising from the interaction of the magnetic field with the spin i particle 

magnetic moment, xb = —fl-B = |— fiB cos tp\ where tp is the tilt angle. For the sake of simplicity, let us set 
[j, = l, then in the transverse case <p = and therefore xb = B = B±. This is our best guess and we justify it by 
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noticing that the magnetic field intensity is indeed a relevant quantity in this experiment (see equation ([73])) and 
its components (intensity) are quantities that are varied during the transitions from integrable to chaotic regimes. 
In the regular regime, we say the magnetic field intensity is set to a well-defined value (xb) = fJ-B- Furthermore, 
notice that the exponential distribution is identified by information theory as the maximum entropy distribution if 
only one piece of information (the expectation value) is known. Finally, the chosen composite probability distribution 
p(mtegrabie) ^ XA ^ XB |^ A; fj,j^j encoding relevant information about the system is given by, 

pCintegrable) ^ ^ = ^(Poisson) ^^(exponential) 



1 

■ exp 



MaMb 







. VMa 


Mb/. 



(77) 



Again, we point out that our probability (|77p is our best guess and, of course, must be consistent with numerical 
simulations and experimental data in order to have some merit. We point out that equation (|77p is not fully justified 
from a theoretical point of view, a situation that occurs due to the lack of a systematic way to select the relevant 
microvariables of the system (and to choose the appropriate information about such microvariables). Let us denote 
^(mtcgiabic) ^ e two-dimensional curved statistical manifold underlying our information geometrodynamics. The line 
element (^ s2 ) intograblc on vVl^ ntcglablc ' ) is given by, 

( ds2 ) integrable = TX^A + "4^1 ■ (78) 

Applying our IGAC to the line element in ([75]) and following the steps provided in the ED Gaussian model of Sections 
II and III of this paper, we obtain polynomial growth in y^ tc s rablc an( j logarithmic IGE growth, 

^ntegrable) (r) ™ ^ ^ntegrable) (r) r^oc + ^ (?g) 

The quantity cjg is a constant proportional to the number of exponential probability distributions in the composite 
distribution used to calculate the IGE and c' IG is a constant that depends on the values assumed by the statistical 
macrovariables [i A and fi B . Equations (|79jl may be interpreted as the information-geometric analogue of the compu- 
tational complexity £)( rcgulai ) ( r ) an d the entanglement entropy <s( rc g ular ) defined in standard quantum information 
theory, respectively. We cannot state they are the same since we are not fully justifying, from a theoretical standpoint, 
our choice of the composite probability ([77]) . 



B. Chaotic Statistical Model: Wigner-Dyson coupled to a Gaussian Bath 

In the chaotic case, the Hamiltonian 7Y( chaotlc ) describes an antiferromagnetic Ising chain immersed in a tilted 
homogeneous magnetic field -Btiitcd = £>_l £>_l + -By -By with the level spacing distribution of its spectrum given by the 
Poisson distribution pO Vl s ner - D y so11 ) (x A \fi' A ) 

(Wigner-Dyson) , / , / >. _ 7TX A I «|\ 

Pa I^aIMa) - 2 ^ 2 cxp ^ ^J. (80) 

where the microvariable x' A is the spacing of the energy levels and the macrovariable /x A is the average spacing. The 
chain is immersed in the tilted magnetic vector field which has two components B± and B\\ in the Hamiltonian 
w (chaotic)_ We translate this piece of information in our IGAC formalism, coupling the probability (|5P[) to a Gaussian 

(Gaussian) / / i / / \ • 1 

Pb (^bIMbj ct b) S lven oy, 

(Gaussian) / I \ I / \ 1 / (^B — Mb) I roi\ 

P ( b >(*b\i*,'*) = ^**<p[ tof-)' ( 81 ) 

where the microvariable x' B is the intensity of the magnetic field, the macrovariable fi' B ^ s the average intensity, and 
a' B is its covariance: during the transition from the integrable to the chaotic regime, the magnetic field intensity 
is being varied (experimentally). It is being tilted and its two components (B± and B») are being varied as well. 
Our best guess based on the experimental mechanism that drives the transitions between the two regimes is that 
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magnetic field intensity ( actually the microvariable fiB cos tp) is Gaussian-distributed (two macrovariables) during 
this change. In the chaotic regime, we say the magnetic field intensity is set to a well-defined value (x' B ) — (a' b with 

covariance ctb = ■J (^(x' B — {x' B )) 2 ^. Furthermore, the Gaussian distribution is identified by information theory as 

the maximum entropy distribution if only the expectation value and the variance are known. Therefore, the chosen 
composite probability distribution p( chaotlc ) ( a; ^ j x' B \[i' A , n' B , a' B ) encoding relevant information about the system is 
given by, 

p(chaotic) / / l\ I I I \ _ (Wigncr-Dyson) ,/ , / v (Gaussian) ,/ , / is 
f (X A , X B \[i A , fl B , <J B ) — p A \ x a\^a)Pb I^bIMb* a B> 



-x A exp 



( x 'b - Mb)' 



2<t' 2 



(82) 



Let us denote _ / v(^ chaotlc ) the three-dimensional curved statistical manifold underlying our information geometrody- 



namics. The line element (ds 2 ) , . onM g 

V / chaotic >-> 



(chaotic) 



is given by, 



( ds2 ) chaotic = TTI^Mchain + ^Mb + ^2<^I 

Ma 



(83) 



Applying our IGAC to the line element in ([83]) and following the steps provided in the ED Gaussian model of Sections 
II and III of this paper, we obtain exponential growth in y^ haotlc ) anc [ li ne ar IGE growth, 



V 



(chaotic) 
M s 



(r) 



C IG exp(JC IG r),S<J ,otic > (r) 



TClGT. 



(84) 



The constant Cjg encodes information about the initial conditions of the statistical macrovariables parametrizing 



elements of M. 



(chaotic) 



. The constant /Cjg, 



IG 



+ 00 dS_\4 s ( T ) 



lim 



-log 

T 



J Ms (0) 



def . 
= A J, 



(85) 



is the model parameter of the chaotic system and depends on the temporal evolution of the statistical macrovariables. 
It plays the role of the standard Lyapunov exponent of a trajectory and it is, in principle, an experimentally observable 
quantity. The quantity Jm s ( t ) 1S the Jacobi field intensity and Xj may be considered the information- geometric 
analogue of the leading Lyapunov exponent in conventional Hamiltonian systems. Given an explicit expression of 
JCig in terms of the observables fj,' A and /i B and a' B , a clear understanding of the relation between the IGE (or 
Kig) and the entanglement entropy (or K, q ) becomes the key point that deserves further study. Equations (|84p are the 
information-geometric analogue of the computational complexity £)( chaotlc ) ( r ) anc j the entanglement entropy ^t 011810 * 10 ) 
defined in standard quantum information theory, respectively. This result is remarkable, but deserves a deeper analysis 
in order to be fully understood. 

One of the major limitations of our findings is the lack of a detailed account of the comparison of the theory 
with experiment. This point will be one of our primary concerns in future works. However, some considerations 
may be carried out at the present stage. The experimental observables in our theoretical models are the statistical 
macrovariables characterizing the composite probability distributions. In the integrable case, where the coupling 
between a Poisson distribution and an exponential one is considered, /i A and fi B are the experimental observables. 
In the chaotic case, where the coupling between a Wigner-Dyson distribution and a Gaussian is considered, n' A 
and n' B and a' B play the role of the experimental observables. We believe one way to test our theory may be 
that of determining a numerical estimate of the leading Lyapunov exponent A max or the Lyapunov spectrum for 
the Hamiltonian systems under investigation directly from experimental data (measurement of a time series) and 
compare it to our theoretical estimate for Aj [65]. However, we are aware that it may be extremely hard to evaluate 
Lyapunov exponents numerically. Otherwise, knowing that the mean values of the positive Lyapunov exponents are 
related to the Kolmogorov-Sinai (KS) dynamical entropy, we suggest to measure the KS entropy K, directly from 
a time signal associated to a suitable combination of our experimental observables and compare it to our indirect 
theoretical estimate for K-ig from the asymptotic behaviors of our statistical macrovariables (66|. We are aware 
that the ground of our discussion is quite qualitative. However, we hope that with additional study, especially in 
clarifying the relation between the IGE and the entanglement entropy, our theoretical characterization presented in 
this paper will find experimental support in the future. Therefore, the statement that our findings may be relevant 
to experiments verifying the existence of chaoticity and related dynamical properties on a macroscopic level in energy 
level statistics in chaotic and regular quantum spin chains is purely a conjecture at this stage. 
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VII. FINAL REMARKS 



In this paper, we reviewed our novel information-geometrodynamical approach to chaos (IGAC) on curved statistical 
manifolds and we emphasized the usefulness of our information-geometrodynamical entropy (IGE) as an indicator 
of chaoticity in a simple application. Furthermore, knowing that integrable and chaotic quantum antiferromagnetic 
Ising chains are characterized by asymptotic logarithmic and linear growths of their operator space entanglement 
entropies, respectively, we applied our IGAC to present an alternative characterization of such systems. Remarkably, 
we have shown that in the former case the IGE exhibits asymptotic logarithmic growth while in the latter case the 
IGE exhibits asymptotic linear growth. 

It is worthwhile emphasizing the following points: the statements that spectral correlations of classically integrable 
systems are well described by Poisson statistics and that quantum spectra of classically chaotic systems are universally 
correlated according to Wigner-Dyson statistics are conjectures, known as the BGS (Bohigas-Giannoni-Schmit, [68| 
and BTG (Berry- Tabor- Gutzwiller, 69]) conjectures, respectively. These two conjectures are very important in the 
study of quantum chaos, however their validity finds some exceptions. Several other cases may be considered. For 
instance, chaotic systems having a spectrum that does not obey a Wigner-Dyson distribution may be considered. A 
chaotic system can also have a spectrum following a Poisson, semi- Poisson, or other types of critical statistics [7Cj. 
Moreover, integrable systems having a spectrum that does not obey a Poisson distribution may be considered as well. 
For instance, the Harper model would represent such a situation. Moreover, it is worthwhile pointing out that not 
every chaotic system characterized by entropy-like quantities growing linearly in time has a spectrum described by a 
Wigner-Dyson distribution. Well-known examples presenting such a situation are the cat maps [7l[ and the famous 
kicked rotator [72| where its spectrum follows a Poisson distribution in cylinder representation and a Wigner-Dyson 
in torus representation but the properties of entropy-like quantities are the same in both representations (at least 
classically). All these cases are not discussed in our characterization. 

Therefore, at present stage, because of the above considerations and because of the lack of experimental evidence 
in support of our theoretical construct, we can only conclude that the IGAC might find some potential applications in 
certain regular and chaotic dynamical systems and this remains only a conjecture. However, we hope that our work 
convincingly shows that this information-geometric approach may be considered a serious effort trying to provide a 
unifying criterion of chaos of both classical and quantum varieties, thus deserving further research. 
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